Static QQ Potentials and the Magnetic Component of QCD Plasma near Tc 
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Static quark-anti-quark potential encodes important information on the chromodynamical inter- 
action between color charges, and recent lattice results show its very nontrivial behavior near the 
deconfinement temperature Tc. In this paper we study such potential in the framework of the "mag- 
netic scenario" for the near Tc QCD plasma, and particularly focus on the linear part (as quantified 
by its slope, the tension) in the potential as well as the strong splitting between the free energy 
and internal energy. By using an analytic "ellipsoidal bag" model, we will quantitatively relate the 
free energy tension to the magnetic condensate density and relate the internal energy tension to the 
thermal monopole density. By converting the lattice results for static potential into density for ther- 
mal monopoles we find the density to be very large around Tc and indicate at quantum coherence, 
in good agreement with direct lattice calculation of such density. A few important consequences for 
heavy ion collisions phenomenology will also be discussed. 

PACS numbers: 12.38.Mh, 25.75.-q, 47.75. -(-£ 



I. INTRODUCTION 

The interaction potential between static quark and 
anti-quark pair is a traditional observable to study the 
quark confineinent mechanism in QCD. It was originally 
inferred from heavy meson spectrum and Regge trajec- 
tories, and has then been extensively studied in lattice 
gauge theories, for reviews see e.g. [i|, ^. Its vacuum 
(T = 0) form is well known, usually represented as a 
sum of a Coulomb part V ~ 1/r, dominant at small 
separation between QQ, and a linear part V = ar dom- 
inant at large separation (see the black solid curve in 
FigH]). The latter implies the confinement of quarks and 
has been interpreted in terms of chromo-electric flux tube 
(or "string") formation between well-separated QQ pair. 
The so-called string tension a in the vacuum (T = 0) has 
been consistently determined by different methods to be 
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With current RHIC and future LHC experimental 
programs exploring excited hadronic matter and quark- 
gluon plasma (QGP) at increasing temperature T, it is 
very important to know the finite T form of the static QQ 
potential, which has recently been calculated by means 
of the lattice QCD, see e.g. [3-[5[. At finite tempera- 
ture, there are actually two potentials associated with 
a QQ pair separated by distance r: one is the free en- 
ergy F{T, r) and the other is the internal energy V{T, r), 
with the difference related to the entropy generated in 
the medium by the QQ pair, i.e. 



V{T, r) = F(T, r)-T [dF/dT) = F{T, r)+TS{T, r) (2) 

What is directly evaluated on lattice is the free energy 
F{T, r) from which the corresponding V{T, r) and S{T, r) 
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FIG. 1: The static QQ potential at T ^ Tc (adapted from 3]). 
The blue boxes are for the internal energy V{r) while the red 
diamonds are for free energy F{r), with the green dashed line 
indicating the strong linear rise in V{r) for r G (0.5, l)fm and 
the black solid line showing the vacuum QQ potential. 



can be inferred [2]. While at T = there is no entropy 
and the free and internal energies are identical, splitting 
between the two shall be expected at T > and may 
carry key information about the medium and deconfine- 
ment transition near Tc. 

The lattice results indeed show remarkably different 
potentials F{T,r) and V{T,r) near Tc (see e.g. Fig. 1-4 
in Q and also here Fig [T] adapted from [3|). In particular 
let us emphasize two important points. 
(i) The tensions (slopes of the potentials at r about 
0.3 — Ifm) have very different temperature dependeces: 
while the tension of the free energy ap decreases with T, 
to near zero at Tc (an expected signal of deconfinement), 
the tension of the internal energy ay remains nonzero till 




FIG. 2: Effective string tensions in the free energy crpiT) 
(from [^) and the internal energy cfv{T) (extracted from [3]). 



about T — l.STc, with a peak value at Tc about 5 times 
(!) the vacuum tension a^ac (see Fig 12]). 
(ii) This drastically different behavior persists to very 
large distances, where linear behavior changes to satu- 
rated values. Near Tc the internal energy flattens to huge 
asymptotic value at large r —>■ cx), e.g. V{T, cxo) ~ 4 GeV 
at Tc with the corresponding entropy S(Tc,(X)) « 20 im- 
plying huge number of states involved, ~ exp{20). 
These features indicate strikingly strong interaction be- 
tween the static color charges and the medium near Tc, 
which persists into the deconfined phase. 

Such static QQ potentials at finite T are closely con- 
nected with a number of phenomenological issues. For 
example, the consequence of these features for the sur- 
vival of quarkonium in deconfined plasma is much de- 
bated, e.g. on what/which potential should be used [3- 
0]. If, as suggested in [ifljllllj the internal energy is used, 
J/tp state would exist even in the deconfined plasma in 
1 — 2Tc. Persistence of some baryonic states above Tc 
is also indicated by other observable like the baryonic 
susceptibilities [1^, [ij|. These potentials also imply sig- 
nificant interaction energy in the quark-gluon plasma and 
in the many body context this may lead to a large clas- 
sical plasma parameter F (defined as the ratio of average 
interaction energy to average kinetic energy) : indeed the 
F value in sQGP has been estimated to be above one 
(about 3) and thus in a typical liquid regime (see for ex- 
ample |14| - KL6| ). If so, QGP would be a strongly coupled 
Coulombic liquid, in agreement with the stron g co Uec- 
tive flow observed at RHIC, see more in reviews pTtilSj. 
Apart from QGP phenomenology, it is important to un- 
derstand the microscopic origin of the potentials, espe- 
cially the strong splitting between two potentials and the 
large energy/entropy associated with the static QQ pair 
near Tc- Earlier attempts can be found in e.g. 19-21]. 

In this paper we will specifically focus on the "ten- 
sions" (Tp and cry (as shown in Figl5]) related to the 



linear part of the potentials (while leaving the discus- 
sion of "screening" behavior at very large distances to 
further studies). We will provide an explanation in the 
framework of the magnetic scenario of QCD plasma near 
Tc [H, [H-ii]. In such a scenario, the near Tc QCD 
plasma is strongly influenced by the magnetic compo- 
nent, made of relatively light and abundant chromo- 
magnetic monopoles. Those are quasiparticles above Tc 
which undergo the Bose-Einstein condensation (BEG) 
below Tc, enforcing color confinement (for reviews see 
e-g- U^ [23 )■ Two key points of the present model for 
the potentials are: First, we identify the QQ free energy 
as been probed by an adiabatically "slow separation" pro- 
cess while the internal energy by a "fast separation" pro- 
cess. Second, we further relate the linear part of poten- 
tials with the fiux tube formation, enabled by condensed 
monopoles below Tc while thermal monopoles above Tc, 
between the QQ pair during the separation process |28| , 
and relate the free/internal energy tensions with the con- 
densed/thermal monopoles respectively. These ideas will 
be elaborated more in Section II and III. 

The rest of the paper is structured as follows. In 
Section-IV we will develop an analytic "elliptic flux bag" 
model for a static charge-anti-charge pair by solving the 
Laplace equation for electric field inside it. This al- 
lows to get the potentials correctly interpolating between 
Goulomb at short distance and linear behavior at larger 
distance. The model will then be used in section-V to de- 
termine the free and potential energies and relate the ex- 
tracted <7f{T) and (Jv{T) with the monopole condensate 
and the thermal monopole density, respectively. Finally 
we summarize the results in Section- VI. 



II. FREE V.S. INTERNAL ENERGY AND 
SLOW V.S. FAST SEPARATION 

Let's start by examining the difference between the 
free energy and the internal energy. We already in- 
troduced the effective string tensions upiT) and cry(T) 
as the slopes of linear parts in F{T,r) and V{T,r) 
respectively, and emphasized their quite different T- 
dependencies shown in Figl2] While crp vanishes at T > 
Tc, (Jy survives to at least l.STc- While ap monotonously 
decreases with T , cry peaks at Tc to a maximal value of 5 
times the vacuum string tension a^ac- What is the differ- 
ence in the meaning of F and V , and why do they have 
such different T-dependence? As has been emphasized 
in [lOj, the free and internal energies actually correspond 
to slow and fast (relative) motion of the charges, respec- 
tively. Let us explain this idea in more details. 

Gonsider the "level crossing" phenomena, occurring 
while the separation between charges is changed. Sup- 
pose a pair of static charges (held by external "hands" ) 
are moved apart in thermal medium at certain speed 
V = L. For each fixed L, there are multiple configu- 
rations of the medium populated thermally. When L is 
changed, the energies of these configurations are cross- 



ing each other, and at each level crossing there is cer- 
tain probability to change between the levels, depend- 
ing on the speed of separation v = L. If the motion is 
adiabatically slow, then all the level crossing processes 
happen with probability 1: in thermodynamical context 
this leads to maintained equilibrium and maximal en- 
tropy/heat generation. If however the pair is separated 
very fast, then the level crossing is suppressed and the 
medium is no longer in equilibrium with the pair. The 
amount of entropy generated is less than in the adiabatic 
case. In the extreme case one may expect that the pair, 
if moving on a time scale much much shorter than the 
medium relaxation time scale, decouples from the media 
and produce negligible entropy. It is plausible, therefore, 
to identify the adiabatic limit as probing the free en- 
ergy F(r, L) measured on the lattice with the presence 
of static QQ pair. The "internal energy" V{T^ L), on the 
other hand, is different from F{T, L) by subtracting the 
entropy term and thus can be probed in the extremely 
fast limit in which possible transitions among multiple 
states via level crossing do not occur and no entropy is 
generated. 

We emphasize that such phenomenon in thermal 
medium is a direct analogue of what exists in pure quan- 
tum mechanical context. Perhaps the oldest example is 
the so called Landau- Zener phenomenon [23, [20] of elec- 
tron dynamics during the vibrational motion of two nu- 
clei in a diatomic molecule. Specific electron quantum 
states ilin{L) are defined at fixed L (the separation be- 
tween two nuclei) with energies En{L), and certain levels 
cross each other at specific value of L. The issue is the 
probability of the transition during such crossing of two 
levels. Consider two levels with their energies given ap- 
proximately by Ei{L) sa aiL + Ci and E2{L) w a2L + C2 
near the crossing point. When the two nuclei approach 
the crossing point adiabatically slowly w = L — >■ 0, the 
electrons always change from one state to the other se- 
lecting the lowest state at any L. If the two nuclei have 
fast relative motion then the transition between the two 
levels at crossing point is suppressed. More quantita- 
tively. Landau and Zener showed that the probability to 
remain in the original state (i.e. no transition) is expo- 
nentially small at small velocity v 



p 



exp 



v\ai ~ (T2I 



(3) 



where H12 is the off-diagonal matrix element of a two- 
level model Hamiltonian describing the transition be- 
tween the two levels. 



III. STABLE AND METASTABLE FLUX TUBES 

We now turn to the possible microscopic origin of the 
linear rise in both potentials. Let's start with the "dual 
superconductor" model for QCD confinement in the vac- 
uum, introduced by t'Hooft-Mandelstam [31| and well 
supported by extensive studies in lattice QCD. In this 



model, certain "magnetically charged" condensate (i.e. 
a magnetic superconductor) occupies the vacuum and 
expels the electric fiux between QQ into a stable fiux 
tube by forming magnetic super- current on tube surface, 
known as (dual) Meissner effect. Such fiux tube naturally 
gives rise to a linear potential and the vacuum string ten- 
sion is thus identified with energy per unit length of such 
fiux tube, mathematically described by the well-known 
Abrikosov-Nielsen-Olesen (ANO) solution 1331 (for re- 
views and further references see e.g. [il,y,|33,|33])- What 
happens at finite T then? With increasing T, the free en- 
ergy tension decreases and eventually the linear part in 
QQ free energy disappears at Tc^ signaling the deconfine- 
nicnt transition. Since the flux tube and free energy ten- 
sion is a direct consequence of the magnetic condensate, 
the decrease of ap toward Tc is naturally interpreted as 
the gradual "melting" of the magnetic condensate due to 
thermal excitations: similar phenomena is known for the 
usual superconductor in condensed matter systems. 

Now, where does the linear part in internal energy (and 
the associated tension ay) come from? In particular, 
why does it p ersist even above Tj.? The answer first pro- 
posed in [1^ relates it to the "normal" monopoles, as 
opposed to the BEC condensed monopoles existing only 
below Tc- Such thermal monopoles can also expel the 
electric fiux into (meta-stable) flux tube by forming a 
magnetic current (which may suffer from dissipation) on 
the tube surface: its dual phenomenon, i.e. magnetic 
flux tube formation in thermal electron plasma is well 
know in classical (e.g. solar) plasma physics. Specific 
condition for the persistence of the electric fiux tube in 
a magnetic plasma was further developed in [23], for in- 
finitely long flux tubes. There it has been found that 
"normal" monopoles are much less effective for this task 
as compared with "super" monopoles, but nevertheless 
able to mechanically stabilize the flux tube provided high 
enough density of these thermal monopoles. What has 
not been previously considered is the mechanism for dy- 
namical formation of flux tube between a QQ pair with 
finite separation. 

Here we provide a dynamical explanation of why large 
energy, growing approximately linearly with length, ap- 
pears in a magnetic plasma when a pair of two electric 
charges are separated with certain speed v, see sketch of 
the setting in FiglJl The answer lies in the Maxwell equa- 
tions (with the presence of magnetic sources, see e.g. 4^]), 
in particularly the dual Faraday's law which relates the 
circulation of the magnetic field / Bdl over a closed con- 
tour with the change of electric flux penetrating the en- 
closed area. As an electric charge moves through the 
loop, rotating magnetic field in the magnetic medium 
leads to solenoidal magnetic current (a "magnetic coil"). 
In the confined phase T < Tc this current, after relax- 
ation, becomes the persistent super-current, remaining 
forever without loss: thus the free energy F has a linear 
term for T < Tc- In a deconfined plasma phase T > T^. 
this is impossible, thus ap = 0: the solenoidal "magnetic 
coil" created in the fast process has only normal magnetic 




FIG. 3: Schematic demonstration of magnetic solenoidal by 
Dural Faraday's law, see text. 



current, which is a meta-stable flux tube and eventuahy 
disappears due to dissipation. Yet it is stih generated: 
thus cry ~ '^F is nonzero and there is sphtting between 
free and internal energy. 

Let us emphasize again the different roles of the super 
and normal magnetic components. The former responds 
quantum mechanically as a whole and does not gener- 
ate any entropy nor contribute to the splitting, the lat- 
ter, however, has finite relaxation time and nonzero dis- 
sipation, "feels" the different time scales involved in the 
slow/fast processes, and therefore is responsible for en- 
tropy generation and the splitting between free and inter- 
nal energy. In short, the ap tells us about the super com- 
ponent only, while the difference ay — crp tells us about 
the normal component. One arrives at the following pic- 
ture for an evolving magnetic medium: with increasing 
T the monopole ensemble starts as a monopole conden- 
sate and continuously evaporates into a mixture of both 
condensed and thermal monopoles; at T > Tc the con- 
densate melts entirely into a normal component of ther- 
mally excited monopoles. If so, the thermal monopoles 
are expected to be most important in the temperature 
range 0.8 — 1.3 Tc where the splitting is most significant. 

We end this Section with discussions on a few impor- 
tant phenomenological implications of the "dual Fara- 
day effect" and the meta-stable flux tube. First, it 
means that magnetic monopoles may induce new mech- 
anism of (electric) jet energy loss particularly near T^. 



In the jet quenching process, a very fast electric parton 
(quark or gluon) penetrates the bulk medium through 
various phases and thus may create behind it the above 
discussed "magnetic coil" in the near Tc region where 
there are abundant thermal monopoles to be accelerated 
solenoidally by B field due to the fast moving electric jet 
and thus take enormous amount of energy away from the 
jet. Such near- Tc enhancement of jet quenching has been 
first suggested by us in [35[ and found to be strongly fa- 
vored by the azimuthal anisotropy data of jet quenching. 
Second, it also implies specific patterns of multi-particle 
correlations if such flux tubes can be created and pro- 
tected by monopoles in heavy ion collisions, as elaborated 
first in ^36i] . One example is related to what happens to 
the flux tube created by a fast jet: clearly the monopoles 
forming the "coil" will subsequently collide with the bulk 
thermal matter, with their energy being converted and 
distributed into the bigger volume: this may possibly 
be the beginning of "conical flow" process suggested in 
[331 • The other example concerns the narrow-azimuthal- 
angle long-range-rapidity correlations known as "ridge" 
which seems originating from certain local initial fluc- 
tuation seeds, but its narrowness in angle may be pos- 
sibly preserved till the end of long bulk evolution only 
if certain mechanism like the flux tube by the thermal 
monopoles protects the initial seed from acoustic expan- 
sion (see detailed discussions in [36J). The existence of 
meta-stable flux tubes in the near T^ plasma (and their 
associated large entropy) may also bear relevance to the 
observed cluster correlations [Sl]. While the existence 
an d dy namical formation of such flux tubes are studied 
in [28] and here, another very important question (par- 
ticularly for phenomenology) is its life time, i.e. the flux 
tube decay and the end products. This problem has re- 
cently been partially addressed in J39J| where a relatively 
short life time is found in the classical treatment. On 
one hand from hydrodynamic modeling we know the near 
Tc plasma has very small shear viscosity which indicates 
short mean free path and frequent scattering, while on 
the other hand for the magnetic currents to last long 
and hold the fiux tube they better do not scatter too 
often: such a dilemma might be resolved if the thermal 
monopoles become really coherent over large distance at 
T close to Tc and the lattice study supports such coher- 
ence [40|. These questions will be studied further else- 
where. 



IV. ELECTRIC FIELD SOLUTION IN THE 
ELLIPSOIDAL BAG 



In this Section we will solve the Maxwell equation for 
electric field induced by a pair of static charge-anti-charge 
separated by a distance L = 2a along z axis [±Qe sitting 
at ^az)^ with a special "tangent boundary condition" 



(T.B.C.) on the boundary surface S^ , i.e. 



suming $ = $(^,77) independent of angle 



^^<t>{r)=Qe[S^{r-az)-S^{r + az)] (4) 

V^-nssIss =0 

The model itself is a version of an old idea known as 
the Bag Model used for hght hadrons [4l| at T = 0, 
now generalized to give an approximate description of 
the electric field configuration between static QQ in the 
chromo-magnetic medium at finite temperature. 

A simplification we use is that the boundary Ss is 
approximated by a rotational ellipsoid with the two 
charges at its focal points. This boundary shape can 
be specified by a single parameter S,b-, the elliptic- 
ity. Such boundary S^ is very conveniently param- 
eterized in terms of the parabolic coordinates system 
(^,77, (/)), which we use: see Appendix A for necessary 
formulae related to it. In Fig|4] we show a few el- 
lipsoidal shapes with parameters (from inside to out- 
side) (L,Cs) to be (0.1, 6.62), (1,1.68), (2, 1.29), (3, 1.16) 
respectively, the dashed lines indicate constant-?] curves 
(for L = 3 case) with (from top to bottom) tj = 
0.8,0.5,0.2,-0.2,-0.5,-0.8, the solid/empty circles in- 
dicate the positions of positive/negative charges, and the 
arrows indicate the tangent electric fields on the bound- 
ary. 
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FIG. 4: The ellipsoidal shapes we use for solving the electric 
field equations, see text for detailed explanations. 



We follow the standard method in classical electrostat- 
ics, see for example [42|. First, we rewrite (|3|) in (^, 77, </>) 
coordinates for solutions with axial symmetry, i.e. as- 
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(5) 



The last line in the above is an expansion of the 77- 
dependence in terms of Legendre functions P^, [77] which 
in the interval 77 G [—1,1] form a set of orthogonal and 
complete basis functions. Similarly, we do the expansion 
for the 77-dependence of $: 



!^=l,3,5, 



ttL 



(2i/ + l)P,[77] 



(6) 



Then by simply comparing the coefficients of Pi^It]] on 
both sides of Eq. ([S]) we obtain the equations for the func- 
tions f^[£_] defined in ^ S (1,00): 



d 
d^ 



t2\'^f'^ 



(l-^^)f 



+ v{v + l)f, = -5{^-\) 



(7) 



while the boundary condition in EqQ now becomes 

/'K = eB] = o (8) 

with the parameter ^b specifying the boundary surface 
Y,B- The solutions are given in terms of the Legendre 
functions of the first and second kinds: 

fAi]^-k^PAC\-QAi] (9) 

P'u[^b] iBPu[^B]-P.-l[^B] 

The full electrostatic potential is then given by 
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(10) 



1^=1,3,5, 



We've used the Neumann expansion of Legendre func- 
tions (see e.g. [43) to write down the second equality: 
in there the first two terms are nothing but the usual 
Coulomb potentials by the iQe charges, while the last 
summation term reflects the nontrivial boundary contri- 
bution. At very large v the summand terms go asymp- 
totically like vC" I^^B^"^ , so with S^ satisfying I < C < ?s, 
the summation is guaranteed to converge. The electric 
field E = — V^ has been calculated using (jA4|) and the 
expression is quite lengthy which we skip showing here. 



The volume occupied by the electric field (i.e. the el- 
lipsoid bulk within ^g) is given by 

Ve{L,Cb) = [ ^ d^ [ dr, f d^H^H^H^ 
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(11) 



And the total electric field energy in this volume is given 
by 
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The Eseif is the familiar self-interaction of the two 
charges which we discard. The "real" interactional en- 
ergy Ee consists (again) a Coulomb piece and a boundary 
modification. 

We conclude this section by one remark: so far the two 
key variables L and fs remain free parameters: they will 
be related in the next section. 



V. THE FREE AND INTERNAL ENERGY OF 
THE CHARGE PAIR 

With the solutions of electric field in the ellipsoidal bag 
(characterized by two parameters L and £,b ) from preced- 
ing Section, we now examine the dynamic formation of 
such bag when separating a pair of QQ from zero to a 
finite distance L. The key point is that for a given L, the 
bag boundary S^b shall be optimized so that the "cost" 
for creating such a configuration is minimized. Further- 
more we study two settings: slow and fast separation of 
the QQ to a finite distance L, with the outcome being 
respectively the free and internal energy associated with 
the pair. With slow separation, the free energy increase 
associated with the pair shall be minimized, and the dom- 
inant contributions to the free energy include both the 
electric field energy stored inside the bag and the energy 
needed to exclude the monopole condensate out of the 
bag volume (noting that for both there is no entropy as- 
sociated and for thermal monopoles their contribution to 
free energy in the slow separation process largely cancels 
out between energy and entropy). With fast separation, 
the energy increase shall be minimized, and the dominant 
contributions to the energy include both the electric field 
energy stored inside the bag and the energy deposited to 



the thermal monopoles via the dual Faraday effect. We 
will calculate both processes in the rest of this Section 
and make connections with the lattice data. 



A. Free Energy from Slow Separation 

As afore-discussed, when the QQ pair is separated in 
an adiabatically slow way, the super component of the 
magnetic medium i.e. the monopole condensate will be 
expelled entirely (in an idealized picture) out of the vol- 
ume Ve occupied by electric field. Suppose the conden- 
sate has a negative energy density —ec (thus a positive 
"bag pressure"), then the overall change in free energy 
brought about by separating the pair will be 



AF =8E{L,^B)+ec[T)-VE{L,^B) 



(13) 



Now for given separation distance L and bulk temper- 
ature T, we determine the physical boundary of fiux 
bag ^3 ^ by minimizing the above Ai^, i.e. the physi- 
cal boundary S,e^{L,T) satisfies: 
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Combining the above with Eq. (|llll2p we then obtain 
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3^1 - 1 rf6 



where we have introduced a length scale 

Ic = {^OiEJl^f-cf'^ 



(15) 



(16) 



with ue = Ql/^TT. This equation could be solved easily 
by numerics. For each L with the above determined 
f^ ^, we obtain via (IT^ the free energy associated with 
the pair as a function of separation L, shown in FiglSja). 
It turns out to be Coulomb at short distance(see the 
magenta dashed curve) plus linear at large distance(see 
the blue dashed line). The occurrence of a linear part is 
due to the physical effect that for large L the medium 
pressure (with which the electric field has to balance) 
limits the transverse size of flux bag (where the field 
gets weak as L increases) to saturate rather than grow 
forever: thus the bag shape approaches a cylinder. 
Mathematically, as L — >■ oo one finds ^^'^ —> 1 but 

L ■ y (fs'^)^ — 1 — > finite. In Fig|3]the four bag shapes 
are at growing L with S,b determined as in the above, 
which clearly shows the shape becomes more and more 
cylindrical at large L. 

By fitting the dimensionless slop of the linear part in 
FiglSfa) we obtain the free energy string tension ap: 



ap = 2.32 X a 



'" X 6i/^ 



(17) 
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FIG. 5: (a)(left) free energy F (in unit of ue/Ic) versus separation L/lc; (b)(Right) monopole condensate energy density 
(ec)^ in unit of yja^ac with the two curves for as being 0.5(upper, red) and l(lower, blue) in Ea. (fT7|) respectively. 



Inversely, since we know ap(T) from lattice as shown 
in Figl21 from the above formula we can infer the T- 
dependence of the monopole condensate energy density 
cc'- see FiglSjb). The two curves are for ap being 0.5(up- 
per) and 1 (lower) respectively. In both cases, ec de- 
creases with T and drops abruptly close to Tc- The inter- 
pretation is natural: toward Tc the monopole condensate 
becomes less and less due to increasing thermal excita- 
tions and eventually dies out around Tc. 

A connection can be made between our result (J17p and 
the dual superconductor model (also known as Abelian 
Higgs model) of vacuum confinement [33J . In that model, 
a quadratic Higgs potential leads to a Higgs condensate 
(the prototype of postulated monopole condensate) (/)o 
(with dimension of mass). By solving ANO flux tube a 
string tension is obtained in the form ^ya = Ci0o with 
the coefficient determined by gauge and Higgs coupling 
constants A and g. On the other hand the Higgs potential 
implies that the condensate has a negative energy density 
—ec — — A(/|q/2, thus one arrives at a similar relation 
between string tension and condensate energy density: 



1/4 



in that model with the coefficient to be 



y/a = C2e^ 

determined numerically for given coupling parameters, 
see e.g. [43]. While that model works primarily at T = 0, 
our model for ap extends to finite T. 



B. Internal Energy from Fast Separation 

Now we study the case of separating the two charges to 
a finite distance L within a time much smaller than the 
relaxation time of the surrounding thermal bath. In par- 
ticular we focus on the region about 0.8 — 1.3Tc, in which 
the normal component of thermal monopoles becomes 
substantial and dominant while the super component be- 
comes less and less. 

During such fast process, each monopole originally 



in the volume to be occupied by the electric filed (i.e. 
the ellipsoidal bag) will get a "kick" due to the dual 
Faraday effect (see Figl3]) but have no time to release 
this energy into the surrounding medium. Suppose the 
positive charge is moved along z-axis from z = C to 
z = C + Sz in 6t (and correspondingly the negative one 
from z = — C to z = — (C + Sz)), then the electric flux 
penetrating the plane z = C changes from to Qg, thus 
generating a magnetic dynamical voltage Qd 5t. For a 
monopole at a transverse distance p from z axis, the 
force is Qm{Q e/ ^t) / {211 p) , thus it gets the "kick" and 
obtains a momentum 5p = QmQe/i^Trp), forming strong 
non-thermal and non-swper magnetic currents. For a bag 
{L,S,b) formed after separation, the total kinetic energy 
passed to the monopoles in the flux bag is obtained by in- 
tegration over the bag volume (with D = QmQ e/4:Tr = 1): 



AK 



M 



4nM D 



L^{e-l){l-ri^) 



n^Dn 



M r2 



l^UCb-i) 



1/2 



(18) 



We emphasize in the above only the monopole density 
nM{T) enters as a property of the medium depending 
on T, while other equilibrium properties of the medium 
shall not be "felt" in such fast process. 

Now the total energy change during such process in- 
cludes the electric fleld energy in the bag volume, the 
energy for expelling the monopole condensate (if there 
is any) out of the volume, and the kinetic energy deliv- 
ered to the normal monopoles, which are summed to be 
AE = £E{L,iB) + ^cVe + AKMiL,(B)- The last new 
term due to the thermal monopoles will then partly con- 
vert into entropy after they interact with the medium 
particles at large for a time longer than the relaxation. 
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FIG. 6: (a)(left) internal energy V (in unit of as/lM) versus separation L/Im\ (b)(Right) thermal monopole density hm /T'^ 
with the two curves connecting box symbols for ue being 0.5(upper, red) and l(lower, blue) in Ea. (|21|l respectively, while the 
two green curves connecting diamond symbols showing S'(7(2)(lower) lattice data and their SC/(3) extrapolation(upper) for 
T > l.STc from [23| (see text for more details). 



ultimately causing the splitting between free/internal en- 
ergy. Since the condensate term ccVe becomes very 
small close to Tc (as we showed in previous subsection) 
and vanishes above T^ we neglect it here for simplicity, 
i.e. AE p:^ £e{L,^b) + AKm{L,^b)- To obtain the phys- 
ical value ^^ ^ , we need to minimize AE according to £,b , 
which leads to 



dAE 



d^E 



This can be written as 



= 



V^ 



IdEi 



2ei - 1 dii 



L 

Im 



(19) 



(20) 



Here we introduced a different length scale Im = 
{2 aE / T^"^ DumY^^ ■ Not surprisingly we find the internal 
energy, shown in FigEJa), to be a Coulomb at short dis- 
tance(see the magenta dashed curve) plus linear at large 
distance(see the blue dashed line). 

Now the string tension in the internal energy is given 
by the following formula: 



1/6 1/3 

X a^ X n^i 



(21) 



Since we know ay from lattice data in the 0.8 — l.STc 
region, by the above formula we can convert aV into 
thermal monopole density nM{T) in the same region: 
see the two curves (connecting box symbols) for ue be- 
ing 0.5(upper, red) and l(lower, blue) respectively in 
FigjSfb). We also show an independent information on 
the thermal monopole density above 1.3Tc from lattice 
study in [25| in FiglBJ^b) as green curves (connecting di- 
amond symbols): the lower one is the original data for 



SU{2) in [25|, while the upper one is an extrapolation to 



SU{2i) by the simple N^—l scaling for monopole species 



(i.e. twice more monopoles in SU{?>) than in S'C/(2)), 
with both curves extended toward Tc according to the 
fitting formula Um/T^ = 0.48/;og(2.48 -T/Tc)^-^^ (and 
twice in the upper one for SU{3)) in [25|. The com- 
parison shows reasonably good agreement between our 
estimates for the thermal monopole density from string 
tension and the measured density by directly identifying 
the thermal monopoles on the lattice. 

A few comments are in order: (i) for 0.8 — 1 Tc the 
density quickly grows toward Tc while at the same time 
results from previous subsection show rapid dropping of 
condensate density in the same region, which strongly 
indicates the scenario that close to Tc monopole conden- 
sate is continuously and substantially getting excited into 
thermal monopoles; (ii) around l.STc we see our results 
connect well to the higher T lattice data with reasonable 
values of coupling a^; (iii) cooling down to Tc we find 
the monopole density um/T^ quickly rising almost by an 
order of magnitude; (iv) the strongly increasing density 
also suggests rapid increase of magnetic screening toward 
Tc, which is in agreement with lattice results [4J|. 

A particularly interesting feature is that the (normal- 
ized) density hm/T^ increases roughly by one order of 
magnitude from l.STc down to Tc, with the number 
much larger than even a Stefan-Boltzman gas. This in- 
dicates that near Tc the monopoles should be very light, 
and furthermore their interactions should make it bene- 
ficial in energy to have a large number of monopole-anti- 
monopole pairs. The monopoles are so dense and light 
that they become the dominant component in the near- 
Tc plasma and presumably become quantum coherent, 
and eventually reach the condensation point at Tc (see 
most recent lattice results in [40| showing evidences for 
such a scenario). It has been suggested that these ther- 
mal monopoles near Tc seem to form a densely packed 
liquid dlltll- 



VI. SUMMARY 

In this paper, we have argued that the free energy 
F{r,T) and internal energy V{rTT) can be probed by 
slow and fast separation of the QQ pair, respectively. 
Furthermore we have identified the linear part in both 
potentials with flux tube formation between the pair: for 
free energy as probed by slow separation, there is sta- 
ble flux tube protected by magnetic super-current due to 
condensed monopoles which has no dissipation and exists 
below Tc] for internal energy as probed by fast separa- 
tion, there is meta-stable flux tube protected by magnetic 
normal current due to thermal monopoles which are very 
dense in the region 0.8 — 1.3 Tc and generate large en- 
tropy (the splitting between free and internal energy) via 
dissipation on longer time scale. Based on these ideas we 
have solved analytically the elliptic bags and provided 
expressions for the potentials at all separations, which 
happen to describe the data very well. 

The main outcome from our study of the static QQ po- 
tentials is the particular relations we suggest between the 
free energy F{r, T) and internal energy V{r, T) measured 
on the lattice and the densities of the condensed and 
"normal" monopoles: see Eqs. lfTTl) and (|21|) . Since those 
densities can be directly obtained from the lattice config- 
urations, one may check if these relations are correct or 
not. Such further tests of the "magnetic scenario" [isl l23l| 
for the near Tc QCD plasma are rather straightforward 
and should be performed. 
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Appendix A 

In this Appendix we briefly list the parabolic coordi- 
nates formulae needed for the calculation in Sec. 2. 

The coordinates we use are (^, 77, </>) with two focal 
points at ±az, which are related to cylindrical coordi- 
nates (p, 4>, z) by 



= a^{e-l){l~V^) , 



<V (Al) 



The variables are defined in the following domains: S, G 

(1,00), 7? G [-1, 1], G [0,27r). Writing ds"^ = H|de^ + 
H'^drj^ + H'^d(tP', we have 



Hf = a 



Ve 



ve 



Hrj = a 






The Laplacian is given by 



(A2) 



-.2 

V = 



a^e-v') 



d_ 

d_ 
drj 
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_^2 _ 1 1 _ jf 

Finally the gradient is given by 



TH^ (A3) 



v = e 



d 



d 



d 



V 



H^dS, Hr^drj ^ H^dcj) 



(A4) 



For more details one could consult books such as HI 
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